Holographic Viscosity of Fundamental Matter 
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A holographic dual of a finite-temperature SU(N C ) gauge theory with a small number of flavours 
Nf <C N c typically contains D-branes in a black hole background. By considering the backreaction 
of the branes, we demonstrate that, to leading order in Nf/N c , the viscosity to entropy ratio in these 
theories saturates the conjectured universal bound r]/s > 1/4-ir. The contribution of the fundamental 



matter 77f und is therefore enhanced at strong 't Hooft coupling A; for example, n fl 



\N c NfT 



four dimensions. Other transport coefficients are analogously enhanced. These results hold with or 
without a baryon number chemical potential. 



PACS numbers: 

Introduction: A universal bound was recently proposed 
1] for the ratio of the shear viscosity to the entropy den- 
sity of any physical system as 77/s > In particu- 
lar, this bound is conjectured to hold for all relativistic 
quantum field theories at finite temperature that exhibit 
hydrodynamic behaviour at long wave-lengths. Perhaps 
surprisingly, experimental results from the Relativistic 
Heavy Ion Collider (RHIC) suggest that, for QCD just 
above the deconfinement phase transition, the value of 
77/s is close to saturating this bound Q- This would in- 
dicate that the quark-gluon plasma formed at RHIC is 
an almost perfect liquid. 

Unfortunately, at present there are no theoretical tools 
with which to calculate transport coefficients in QCD in 
this strong coupling regime, e.g., the viscosity. However, 
a large class of gauge theories are accessible to study with 
the gauge/gravity correspondence [J. In particular, in 
the gauge theory limit of large N c and large 't Hooft cou- 
pling A, the dual description reduces to classical super- 
gravity. In fact, the proposal Q for a universal bound 
on 77/s originated with calculations in this holographic 
context. Explicit calculations 0,|EE10 anc ^ general ar- 
guments 0, M El have demonstrated that the bound is 
exactly saturated by a large class of holographic theories 
in the limit cited above. In order to make contact with 
real- world QCD, it is clearly important to consider 1/A 
and l/N c corrections. For four-dimensional Af = 4 su- 
per Yang- Mills (SYM), the leading correction of the first 
typ e was shown to raise the value of 77/s above the bound 
jlfjj . That is, at finite coupling the bound still holds but 
is no longer saturated. 

A feature common to all of the gauge theories consid- 
ered in these hydrodynamic studies is that the matter 
degrees of freedom transform in the adjoint representa- 
tion of the gauge group 0]. In this letter, we study 
the effect of adding matter fields transforming in the 
fundamental representation, bringing us one step closer 
to QCD. In particular, we focus on four-dimensional 
SU(N C ) SYM coupled to N { fundamental hypermulti- 
plets with N f <C N c . Large- N c counting rules imply that, 
in the deconfined phase, the contribution of the gluons 



and adjoint matter to physical quantities is of order A 2 . 
Further, the first correction in the absence of fundamen- 
tal matter is of order 1, i.e., the relative contribution 
is suppressed by 1/iV 2 . Instead, the relative contribu- 
tion of fundamental matter is only suppressed by N t /N B , 
and therefore it constitutes the leading correction in the 
large- N c limit. 

The dual gravity description is given by N t D7-brane 
probes 0] m the background geometry of N c D3-branes. 
At finite temperature, the latter contains a black hole 
[l4|. We will show that, to leading order in N f /N a , the 
calculation of the 77/s ratio can be effectively reduced 
to one in five-dimensional Einstein gravity coupled to a 
scalar field. General results 0, @ then guarantee that 
77/s = 1/4-7T. Since the D7-brane contribution to the en- 
tropy density is known to be of order s fund ~ XN c N t T 
|15| . this implies that the contribution of the fundamen- 
tal matter to the shear viscosity at strong 't Hooft cou- 
pling is enhanced with respect to that dictated solely by 
large- N c counting rules. 

In the final section we will argue that an analogous 
enhancement takes place for other transport coefficients. 
We will also explain how our results extend straightfor- 
wardly to other holographic gauge theories described by 
Dp/Dq [H |3 or Dp/Dq/Dq O systems, as well as 
to systems with a non-zero baryon number chemical po- 
tential. We will finish with some comments on effects 
beyond order N f /N c . 

Holographic Framework: The shear viscosity of the 
gauge theory in a two-plane labelled by x l ', x 3 may be ex- 
tracted from the retarded correlator of two stress energy 
tensors via Kubo's formula 
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([7y(z),7y(0)]) , (1) 



where no summation over i,j is implied. The stress en- 
ergy tensor is dual on the string side to a metric per- 
turbation Hij polarised along the same two-plane. The 
two-point function above may be calculated by taking 
two functional derivatives of the on-shcll string effective 
action with respect to this perturbation 0] . In the large- 
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N B , large-A limit, this effective action reduces to the type 
1IB supergravity action coupled to the worldvolume ac- 
tion of the D7-branes, I = I UB + J D7 . Schematically, we 
have: 



1 



16ttG 



d w x\f--gR — N { T D7 I d s xy^g~ d + 



(2) 

where g lnA is the induced metric on the D7-branes. In 
terms of the string length and coupling: 



16ttG = (2tt) 't s gi , T D7 = 2tt / \2nt a f g s 



(3) 



The ratio between the normalisations of the two terms 
above is 



e = 16^GA f T D7 = A A 
Ztt iV c 



(4) 



where A = g 2 K N c = 2irg s N c is the 't Hooft coupling. This 
ratio controls the relative magnitude of the D7-branes' 
contribution to physical quantities, e.g., the entropy den- 
sity 0|- We will assume that e<l and hence that the 
D7-branes can be treated as a small perturbation; for 
fixed A this is achieved by taking N t <§; N c . We will 
begin by examining contributions of order e in the next 
section. In the last section, we will comment on effects 
of order e 2 and higher. 

In the absence of D7-branes, the supergravity back- 
ground dual to four-dimensional AT = 4 SYM at temper- 
ature T is (in the notation of |l5j l 



ds 2 
ds 2 



ds 2 + L 2 dn 2 5 , 
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where 



f(p) = 1-1/ p\ f{p) = l + l/p i 



(5) 
(6) 

(7) 



and L = (47rg s iV c ) 1 / 4 ^ ;s is the asymptotic AdS radius. 
There are also N c units of Ramond-Ramond (RR) flux 
through the five-sphere while the remaining supergravity 
fields vanish. The metric ijoKill possesses an event horizon 
at p = 1 . The entropy density of the gauge theory is then 
given by the geometric entropy of the horizon 
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N 2 T 3 , 



(8) 



where G 5 = G/tt 3 L 5 is the five-dimensional Newton's 
constant obtained by dimensional reduction on the five- 
sphere. 

Now we introduce the D7-branes oriented such that 
five worldvolume directions match those of the five- 
dimensional black hole ©, y = {t, x l ,p}, and the re- 
mainder wrap an S 3 (with a possibly varying radius) in- 
side the S 5 of (JSJ. We adapt coordinates in this internal 
space such that 



dni = d6 2 



sin 2 6dnj 



(9) 



and describe the D7-branes embedding as x = x{p)-> with 
X = cosO. To order e°, this is determined by extremising 
the D7-brane action in the background JSJSJ. Asymptot- 
ically, one finds 



X 



m 
P 



c 



(10) 



where m and c are proportional to the quark mass M q 
and condensate (tpip)y respectively. In the interior, the 
D7-branes may or may not reach the black hole horizon 

m 

Viscous Branes: As alluded to above, the calculation of 
the shear viscosity proceeds as follows. First, one solves 
the (linearised) equation of motion for a metric pertur- 
bation H around the appropriate background. Next, one 
evaluates the appropriate action for the perturbed back- 
ground to quadratic order in H. A second derivative of 
the on-shell action then yields the desired two-point func- 
tion with which the Kubo formula (JJJ is evaluated. 

In the absence of D7-branes the appropriate action is 
J IIB and the background is given by JSJ- In the presence 
of the D7-branes the relevant action is instead that in 
cq. J2J). To first order in the £-expansion, this affects 
the calculation of the viscosity in three ways. First, the 
branes will produce 0(e) corrections to the metric (0, 
as well as to the dilaton and the RR axion, since they 
act as new sources in the field equations arising from 
the combined action. These background corrections then 
lead to modifications of the field equation satisfied by 
H. Second, the branes will also modify the H field equa- 
tion directly through the extra source terms originating 
from the variation of I U7 with respect to H . Third, the 
second-derivative of the on-shell action, which yields the 
correlator in Kubo's formula, may acquire contributions 
from/ D7 . |3 

We will now show that only the first two types of pos- 
sible modifications do actually contribute and, moreover, 
that the only type of background correction that needs 
to be considered is the zero-mode (on the five-sphere) of 
the five-dimensional black hole metric ©• 

We begin by considering the third set of possible contri- 
butions listed above. Expanding the brane action around 
the O(e ) background, one finds that, because H en- 
ters the action non-derivatively, the H 2 terms do not 
have a form which will contribute in the Kubo formula 
[Isj . However, turning on H also induces a correction 
Sx = 0(H 2 ) in the embedding of the branes. This leads 
to a surface term in the variation of the D7-brane action, 



dC D7 
d(d pX ) 



(11) 



of the right form to contribute to the two-point corre- 
lator. However, this contribution is proportional to 5m 
[l5j and hence vanishes because the variation of the ac- 
tion with respect to H must be taken while keeping the 
quark mass fixed. 
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Consider now corrections to the background J5J . While 
there are 0(e) corrections to the dilaton and the RR 
axion, these only produce 0(e 2 ) contributions to the H 
field equation because they enter the supergravity action 
quadratically (in the Einstein frame). We are thus left to 
consider the contributions of corrections to the spacetime 
metric. Schematically, to order e, the background metric 
takes the form 



9 = 5o (y) +£9i\y) 



(y)YW((l), (12) 



where y«(fl) are spherical harmonics on S 5 . The 
zero'th order metric 10 and the correction are con- 
stant on the S 5 . However, since the D7-branes only fill 
an S* 3 in the internal space, they also source the g± cor- 
rections which vary over the S 5 . For the following, an 
important point is that the functions gf (y) and g± (y) 
respect the symmetries of the background geometry l|Sjl 
and the brane embedding, i.e., translations in {t, x 1 } and 
50(3) rotations in x l , as well as SO(4) rotations in the 
internal S 3 wrapped by the D7-brane. The perturbation 
H has a similar decomposition: 



H = H {y) 



:J2H[ £ \y)YW(n). (13) 



In the absence of D7-branes, one works consistently with 
the S^-independent perturbation H 0,0- However, in 

the presence of the D7-branes, nontrivial H{ are sourced 
when Hq is turned on. Indeed, after integration over the 
S 5 , the supergravity action produces couplings of the 
schematic form e 2 J d 5 y HqH^ gf . Similarly, the D7- 
branes action produces couplings like e 2 J d b y HqH^ , 
which arise from spherical harmonics that are con- 
stant on the S 3 wrapped by the D7-branes. However, as 
indicated, both types of terms are of order e 2 and so we 
may neglect their contribution since here we only wish to 
determine the correlator (JIJ up to order e. 

We therefore conclude that, to order e, we need only 
consider the zero-modes gf\y) and H^\y), and may 
drop all modes with non-trivial dependence on the S 5 
directions. Hence in working to order e, evaluation of 
the viscosity actually reduces to a five-dimensional cal- 
culation. We can make the latter concrete by dimension- 
ally reducing the action J5| to five dimensions ignoring 
all the Kaluza-Klein modes on the five-sphere, as well as 
the other supergravity fields. By the previous arguments, 
the resulting action still captures all of the relevant fields 
for the calculation of the viscosity to this order. The 
five-dimensional action can be written as 




where g denotes the five-dimensional metric. The first 
two terms originate from the reduction of 7 IIB , whereas 



the last one comes from the reduction of I D7 15]. Note 
that, in the action (|14fl . we have only allowed scalar 
field configurations depending on the radial coordinate, 
since this suffices for our purposes. This system is just 
five-dimensional Einstein gravity coupled to a cosmolog- 
ical constant and a(n unusual) scalar field x- I n an £ ~ 
expansion, the black hole solutions generated by this aux- 
iliary theory will match the asymptotically AdS part of 
the original ten-dimensional solution to order e, i.e., the 
brane profile x{p) an d the background metric l@ plus 
the order-e correction g^ (p) . 

The viscosity may now be obtained by calculating the 
perturbation Hij around the five-dimensional solution 
and taking the second functional derivative of the ac- 
tion (|14fl evaluated on-shcll. However, the black hole 
solutions of our auxiliary five-dimensional system satisfy 
the symmetries required in and hence the result is 
guaranteed to satisfy ij/s = 1/Att. We thus conclude 
that this universal bound is still saturated in the full ten- 
dimensional string theory when working to first order in 
e. An immediate consequence is that the contributions 
of the fundamental matter to the viscosity and the en- 
tropy density are related (within our approximations) as 
^tund = Sfund/47r. The leading contribution to the entropy 
density was determined 15] to be 
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N c NfT h 
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(15) 



where the function h(x) satisfies h(0) — 0, h(oo) — 1, and 
makes a cross-over between both values around x ~ 1. 
Note that this cross-over includes a small discontinuity 
arising from a first-order phase transition of the funda- 
mental matter . We therefore conclude that both s fund 
and 77f und are enhanced at strong 't Hooft coupling with 
respect to the 0(N c N ! )-value dictated solely by large- iV c 
counting rules. 

The calculation of s fund in |15| was performed by identi- 
fying the Euclidean action of the D7-branes with F iund /T, 
where F fund is the free energy contribution of the fun- 
damental matter. The entropy is then determined as 
Sfund = —dF land /dT. This entropy should, of course, co- 
incide with the change in the horizon area induced by 
the presence of the D7-branes. The latter can be ex- 
plicitly verified for the case of massless quarks, which 
corresponds to an 'equatorial embedding' with x — 0. In 
this case, the result from [l5j for the entropy density is 
Sfund = AA c A f T 3 /16. In the action (O, we see that the 
net effect of these 'equatorial' D7-branes is to shift the 
effective cosmological constant. The corresponding black 
hole solution is still given by I©, with the replacement 
L 2 — * L 2 /(l —e/Qnr). The same replacement in (JHJ shifts 
the entropy to order e by 5s = XN C N { T 3 /16, in perfect 
agreement with the previous result. 
Discussion: We have seen that the calculation of the 
contribution of fundamental matter to the shear viscosity 
may be effectively reduced to a calculation in five dimen- 
sions. An analogous simplification takes place for other 
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transport coefficients that can be extracted from corre- 
iators involving local operators with vanishing R-charge, 
since these are dual to modes that carry no angular mo- 
mentum on the S 5 . Examples involving components of 
the stress-energy tensor include the speed of sound v s 
and the bulk viscosity £. Other transport coefficients 
that involve R-charged operators, such as the R-charge 
diffusion constant |5j, or extended strings, such as the jet 
quenching parameter q (see e.g., 0|), may require a ten- 
dimensional calculation. Generically, however, we expect 
the relative contribution of the fundamental matter to be 
of order e ~ XN { /N C> since this controls the backreaction 
of the branes. A special case is the bulk viscosity, which 
may be extracted from a two-point function of T„ — v 2 Tqo 
[22| . This combination vanishes if N f = by conformal 
invariance, and hence £ = 0{e 2 ) in the presence of fun- 
damental matter. 

Above, our discussion focussed on the D3/D7 system, 
but the arguments are easily extended to a more general 
system of Dq-branes in a Dp-background. In particu- 
lar, systems arising from Dp- and Dq-branes intersecting 
over d common spatial directions have been of some in- 
terest 0, |23| . These constructions are dual to a finite- 
temperature SYM theory in p + 1 dimensions coupled 
to fundamental matter confined to a (d+ l)-dimensional 
defect. One new feature in these generalised configura- 
tions is that the defect breaks translational invariance 
along the p — d orthogonal directions. In order to cal- 
culate the shear viscosity [24| along the translationally 
invariant directions parallel to the defect, the simplest 
approach is to compactify these extra directions [26j . The 
arguments in the previous section go through essentially 
unchanged except for the fact that the index £ now labels 
momentum modes both along the S 8 ~ p transverse to the 
Dp-branes and along the p — d directions orthogonal to 
the defect. In this case the problem of calculating the 
leading contribution of the fundamental matter to the 
viscosity /entropy ratio can be reduced to a calculation 
in (d + 2)-dimensional Einstein gravity coupled to a set 
of scalar fields. In addition to the scalar \ above de- 
scribing the size of the internal S n of the Dq-brane, this 
set now includes the dilaton and the metric components 
governing the size of the internal 5 8 ~ p of the background 
geometry and the size of the (p — <i)-dimensional space 
orthogonal to the defect [27j ■ This lower dimensional the- 
ory again captures all of the relevant fields to calculate 
the viscosity to leading order in N f /N . Further, the form 
of the (d + 2)-dimensional gravity theory and the back- 
ground guarantees that rj/s — The leading result 
for the entropy density was determined in [l5| and hence 
we have 

77 fund ~ N N t T d g c!e (T) ^ , (16) 

where g 2 it {T) = \T P ~ 3 is the dimensionless effective 't 
Hooft coupling for a (p + l)-dimensional theory at tem- 
perature T 28]. Here the gauge/gravity duality is only 
valid in the strongly coupled regime [28| and hence we 
see again an enhancement beyond the large- N counting. 



It would be interesting to understand if this enhance- 
ment extends to other transport properties in the same 
way, as was found for the D3/D7 case. The same line 
of argument can also be implemented for the Dp/Dq/Dq 
systems which have also been studied recently 17] . 

The U(N t ) ^ SU(N t ) x U(1) B gauge symmetry on the 
Dq-branes is a global, flavour symmetry of the dual gauge 
theory. The results above also hold when a baryon num- 
ber chemical potential for the U(1) B charge is introduced. 
This is dual |29j] to turning on the time component of the 
gauge potential on the Dq-branes, Ao(p) (30|. The argu- 
ments of the previous section again go through essentially 
unchanged, except for the fact that an additional vector 
is added to the reduced (d + 2)-dimensional Einstein 
gravity theory. Thus the saturation of the bound is not 
affected by the introduction of a chemical potential. 

Above we have worked to the lowest order in the pa- 
rameter e ~ \N;/N C that controls the backreaction of the 
D7-branes on the D3-brane geometry. We have argued 
that to this order one may ignore all effects of this back- 
reaction except for those on the non-compact part of the 
metric. We regard the agreement between the entropy 
density as calculated in | l5| and as obtained here from 
the change in the horizon area as a quantitative consis- 
tency check of this approach. In calculating beyond order 
e, the internal modes, e.g., g[ e \y), and other supergrav- 
ity fields, e.g., the RR axion, will all play a role. Further 
at <D(e 2 ), quantum effects will have to be considered [3l| . 
Closed string loop corrections naturally appear in an ex- 
pansion in g 2 ~ X 2 /N 2 . Thus if as above N t is fixed, 
the loop corrections may be of the same magnitude as 
the higher order D7-brane contributions, i.e., g 2 ~ e 2 if 
N f = 0(l). 

In order to suppress string loop corrections with re- 
spect to the backreaction of the D7-branes (or more gen- 
erally the Dq-branes), one must keep N t /N c fixed while 
taking N c — > oo. In this limit, the entire backreaction 
of the branes must be taken into account. A fully back- 
reacted solution is singular for the D3/D7 system [32| . 
since the dual gauge theory possesses a Landau pole at 
some finite scale A uv . This of course should not affect 
the hydrodynamic behaviour as long as the other scales 
in the problem, the temperature and the quark mass, 
are much lower than A uv (just like the transport prop- 
erties of an electromagnetic plasma are not affected by 
the Landau pole in QED). Similarly, in the case of a gen- 
eral Dp-brane background (with p ^ 3), one finds strong 
coupling or curvature divergences in the far UV, which 
are irrelevant for the long-wavelength hydrodynamics. In 
these cases, the classical backreacted solution would ef- 
fectively resum the effects of the fundamental matter to 
leading order in the 1/A, 1/N C expansions. It may be pos- 
sible to prove that the viscosity bound is still saturated 
by these solutions by extending the ar gum ents of 0, . 
The finite-temperature backgrounds of |L2ll33l ] represent 
a step towards this goal. 
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